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ABSTRACT 

This review describes. Interprets, and compares eleven sensitivity 
experiments that have been made with general circulation models to see how 
land— surface boundary conditions can influence the rainfall, temperature, 
and motion fields of the atmosphere. In one group of experiments, different 
soil moistures or albedos are prescribed as time-invariant boundary conditions. 
In a second group, different soil moistures or different albedos are initially 
prescribed, and the soil moisture (but not the albedo) is allowed to change 
with time according to the governing equations for soil moisture. In a 
third group, the results of constant versus time-dependent soil moistures 
are compared. 

All of the experiments show that the atmosphere is sensitive to the land- 
surface evapotransplratlon: so that changes in the available soil moisture or 
changes in the albedo (which affects the energy available for evapotransplra- 
tlon) produce large changes in the numerically simulated climates. 


*Revlew paper presented at the JSC Study Conference on Land Surface Processes 
in Atmospheric General Circulation Models, Greenbelt, U.S.A., 5-10 jJanuary 1981 . 
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Some Observational and Theoretical Consideration!^ 

Averaged for Che globe and for the year, Che measured river wacer drainage 
from Che condnenca la about a Chlrd as large aa Che aeaaured precipitation 
(Baumgartner and Reichel, 1975, Table 9; Korzun, 1978, Table 150). Thla eeana 
Chat on the average the land-aurface evapotranaplratlon is about two-thirds aa 
large a8 Che preclplcadon. 

In some contlnenCal regions, during pare of Che year, Che evapotransplra- 
cion la larger chan Che precipitation. Thla cannot be known from measurements 
of river flow, but can be derived from measurements of Che CransporC of water 
vapor by Che aCmosphere. An example of Chls for Che cencral and eastern United 
States, in July, is shown in Fig. 1. 

On the left in Che figure is the vertical distribution of the water vapor 
transport divergence, as given by twice daily rawinsonde measurements for two 
July months, and averaged for the region 80*W-100*W, 30*N-47.5*N, which is 
an area of about (2000 km) 2 (Rasmusson, 1968, Table 1 and Fig. 2). From 
Che surface to the 930 mb level there is a water vapor CransporC convergence 
of 1.4 gm cm -2 month -1 : or 14 mm/mo equivalent water depth. Above the 

930 mb level there is divergence of 36 mm/rno. Integrated over Che entire 
depth of the atmosphere there is a net divergence (a net removal of water 
from the region) of 22 mm/mo. 

From the beginning to the end of July the change in the water vapor con- 
tent of the atmosphere is very small. Therefore, the 22 mm of water Chat are 
removed from the region must come from the water stored in the soil; which 
means Chat the evapotransplration is 22 mm/mo larger than the precipitation. 
Inasmuch as the measured average July precipitation in Chls region is about 
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Fig. 1 Water budget (na/no) for the central and eastern United States 
in July. Divergence of the water vapor transport is fron Rassusson (1968). 



94 aa/ao, the average July evapotranaplration must be about 116 tm/mo 
(3.7 aa/day). {A comparable analysis for the central and eastern United 
States by Benton et al. (1953, Figs. 24, 26) gave a July evapotransplratlon 
of 121 aa/ao (3.9 aa/day)]. We can Interpret this water budget as follows: 

The July net radiational heating of the ground in this region, in 
July, is about 140 watt/a 2 (Budyko, 1963, Plate 21), which if used entirely 
for evapotransplratlon would put about 150 aa/ao (4.8 aa/day) of water 
into the air. But if we accept the serologically derived evapotransplratlon 
of 3.7 aa/day (LE - 107 watt/m 2 ), there will be a sensible heat transfer 

froa the ground to the atmosphere, H - (Rjj - LE) - (140 - 107) ■ 33 

and 

watt/a 2 ; .a Bowen ratio, (H/LE), equal to 0.31. 

w\ 

The 116 aa/ao of water vapor, forced into the atmospheric planetary 
boundary layer by the radiational heating of the surface, combines with the 
14 aa/ao brought Into the region by the water vapor transport convergence in 
the boundary layer; and the total of 130 aa/ao of water vapor are transferred 
froa the boundary layer to the free atmosphere. 

In July, the condensation and precipitation In this region is predom- 
inantly of the convective type, with relatively little large-scale upglide 
condensation and precipitation. Therefore, the transfer of water vapor 
upward froa the surface is predominantly by small scale turbulent alxlng 
within the planetary boundary layer, with a handover to cuaulus convection 
which carries the water vapor froa the top of the boundary layer into the 
free atmosphere. 

Of the 130 aa/ao of aolsture carried into the free ataosphere by the 
cuaulus cloud towers, 36 aa/ao (in the fora of water vapor, liquid water 
droplets and ice crystals) are detrained froa the clouds into the cloud 
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environment (where the water droplets and ice crystals evaporate) and are 
removed from the region by the divergence of the water vapor transport in 
the free atmosphere. The remaining 94 ma/mo return to the earth's surface as 
the convective precipitation. The excess of the evapotransplratlon over the 
precipitation , 22 am/mo, Is the moisture withdrawn from the soil. 

According to this analysis, the convective precipitation draws all of Its 
moisture from the water vapor In the planetary boundary layer; and the amount 
of water vapor supplied to the boundary layer by the surface evapotransplratlon 
is an order of magnitude larger than the amount supplied by the water vapor 
transport convergence* This suggests that the surface evapotransplratlon la 
the main determinant of the precipitation. 

[The winter season water budget over the central and eastern United States 
is very different from that shown In Fig. 1. In winter the water vapor 
transport convergence does not change sign with height, but la convergent at all 
levels and produces a net import of water vapor to the region (Rasmusson, 1968, 
Table 1 and Fig. 2). In winter, the condensation and precipitation is pre- 
dominantly of the large-scale upglide condensation type (frontal cloud and 
precipitation) which draws from the water vapor at all levels in the tropos- 
phere. Moreover, In winter the net radlatlonal heating of the ground la 
small (Budyko, 1963, plate IS) and, conaequently, over the unforested 
part of this region the evapotransplratlon la small. In winter, therefore, 
the land-surface evapotransplratlon cannot have much Influence on the 
precipitation or other fields. It Is only In the tropics and In the summer 
season extratropics, where evapotransplratlon Is large and where the 
precipitation Is of the type thst draws Its water vapor from the planetary 
boundary layer, that the land-surface evapotransplratlon can be of major 
Importance. ] 
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With respect to the tropics end the summer season extratropics, two 
questions Immediately come to mind: 

1) If the surface evapotransplratlon, by some means, la greatly 
reduced, can the boundary layer water vapor transport convergence increase 
by a corresponding amount and, in that way, maintain the precipitation? And, 
il) If the surface evapotransplratlon and the boundary layer water vapor 
transpo.t convergence remain the same, can an Increased detralnment and water 
vapor transport divergence in the free atmosphere stop the precipitation? 

e To answer the first question, we write the water vapor transport 
convergence as -V»q v ■ - v*^q — qV* v» where q is the water vapor mixing 
ratio and \v Is the horizontal velocity of the air* 

- v*Vq Is positive when the air that leaves the region is drier than 
the air that enters. But If this drying Is due to the removal of water 
vapor from the boundary layer by cumulus convection, then the convection 
will stop as soon as there Is a small reduction in the boundary layer water 
vapor content. - v*7q cannot maintain the observed rate of precipitation 
over a distance which Is greater than just a few cumulus convection cells; 
say, a total distance of a few kilometers. 

The other term, -qV* \v, also can maintain the observed rate of precipi- 
tation only over a restricted domain in the extratropics. The water vapor 
mixing ratio, in the boundary layer of a maritime tropical air mass over 
the extratroplcal continents in summer, is of the order of 10 parts per 
thousand. Therefore, a boundary layer that is 100 mb (1 km) deep must have a 
horizontal velocity convergence, -V* w, of 0.37 day”* (0.43 x 10"^ sec - *) to 
produce a water vapor transport convergence of 3.7 mm/day. The characteristic 
velocity of the boundary layer air in the extratropics, in summer, is 2 to 
3 n/s; and the angle between this vector velocity and the streamline of the 
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non-dive rgent flow, Integrated over the depth A the boundary layer, la about 10*. 
Thua, if we conalder the circular region with radlua r, we have -?• v ■ 

0.43 x 10“ 5 a*" 1 -(3 b a -1 aln 10*) 2»r/»r 2 ; or r - 240 ka, and 480 ka la 
the Halting dlaaeter of the region where water vapor tranaport convergence 
in the boundary layer can produce the obaerved rate of precipitation. 

In the extratropica, therefore, there la a alte Halt, of the order 
of a few hundred klloaetera, beyond which boundery layer water vapor 
tranaport convergence cannot coapenaate for dlainlahed evapotranaplratlon. 

It la only near the equator, where the divergent coaponent of the velocity 
field la larger and the planetary boundary layer la deeper, that there can 

be appreciable water vapor tranaport convergence over a auch larger elced 
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ee The anawer to the aecond question, "Can an Increaaed detralnaent and 
water vapor tranaport divergence In the free ataoaphere atop the precipita- 
tion?", dcpenda on whether the free ataoaphere la supplied with dry air into 
which the cuaulus cloud towera can detrain. That will happen only If, In 
addition to the boundary layer aasa (and water vapor) convergence, there la 
also a aass convergence In the upperaost troposphere. Then, the cuaulus 
cloud towers can detrain all of the water into the subsiding and diverging 
dry air of the aiddle troposphere, and not produce any precipitation at all. 

The best known exaaple of extensive fields of non-preclpitatlng cuaulus 
clouds of thla kind are the Trade Wind cuaull over the tropical oceans, 
where the subsiding air in the aiddle troposphere has its origin In the high 
level outflow above the lntertroplcal convergence sone. He also see auch 
fair weather cuaulus clouds reaovlng water vapor froa the boundary 
layer, without producing precipitation, west of the trough lines and east 
of the ridge lines of the fast-transient and slow-transient waves in the 
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extratroplcal westerlies, where both the longitudinal and the latitudinal 
acale can be as large aa a few thouaand klloaeters. We can say, therefore, 
that on a scale larger than a few hundred klloaeters jin the extratroplcSj 
land-surface evapotranspl ration is a necessary (but not sufficient) condition 
for convective precipitation. The upper tropospheric circulation aust also 
be favorable for precipitation. 

Because so aany interactive theraodynaalcal and hydrodynaalcal processes 
are Involved, the best way to deteralne the overall Influence of the land- 
surface boundary conditions on the rainfall, teaperature and circulation la 
through experiaents with ataospheric general circulation aodels. Existing 
general circulation aodels have been fairly successful In slaulating the 
observed diaate of the earth, Including the principle geographical and 
seasonal characteristics of the precipitation (WMO , 1979). By Baking pairs 
of tiae-lntegratlons , with all of the Initial conditions and boundary 
conditions the saae except for those which can affect the land-surface evapo- 
transplration, and coaparlng the two solutions, we can ascertain what the 
land-surface Influence Is. 

In the existing general circulation aodels, the two boundary conditions 
that can affect the land-surface evapotranspl rat Ion are the soil aolsture and 
the surface albedo The soil aolsture deteraines how large the evapotranspl ra- 
atlon will be relative to the aodel calculated potential evapotransplratlon 
(the evapotransplratlon when soil aolsture Is fully available): the albedo Is 
a aajor factor In deteralnlng the potential evapotransplratlon Itself. 

The experiaents that are being reviewed are grouped as follows: 

I. Experiaents with non-lnteractlve soil aolsture. 

II. Experiaents with Interactive soil aolsture. 

III. Hybrid experiaents. 
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!n the first group, cither different soil moisture sveilebilities or 
different elbedos are prescribed, and both of these parameters are kept 
constant with time. Such experiments reveal the sensitivity of the atmosphere 
to the boundary conditions. (These experiments are analagous to sensitivity 
experiments in which different non-lnteractlve ocean surface temperatures are 
prescribed: the so-called sea surface temperature anomaly experiments). 

In the second group, the soil moisture (but not the albedo) la inter- 
active and changes with time according to the model's governing equations for 
soil moisture. When the albedos are the same in a pair of comparison runs, 
but the initial soil moistures are different, the Integrations will either 
produce time-series that remain separate (intransitive) or converge to a 
common solution; and, if transitive, they will show how long it takes for 
the two initially different states to converge to a common state. When the 
albedos are different, this trill be another kind of sensitivity experiment. 

In the third group, the hybrid experiments, calculations with non- 
lnteractlve and interactive soil moistures are compared. To the extent 
that the calculation with interactive soil moisture simulates the observed 
rainfall, temperature and circulation of the earth's atmosphere, the com- 
parison will show how the earth's climate may be affected by such imposed 
changes in the land-surface evapotransplratlon as might be brought about by 
large scale deforestation or afforestation, by soil erosion or reclamation, 
or by large scale irrigation. 
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I. Exparlaa nta with Non-Intaractlva Soil Molatura: 

A. Dlfftrtnt soil aolaturaa, with hm albado. 

(1) Shuk*a and Mints (1981). 

(2) Suaraz and Arakawa (paraonal comaunlcntlon) . 

(3) Mlyakoda ami Strlcklar (1981). 

B. Dlffarant albadoa, with «aw soil wolatura. 

(4) Charnay, Quirk, Chow and Kornflald (1977). 

(5) Carson and Sangstar (1981). 

II. Esparlwanta with Intaractlva Soil Molatura: 

A. Dlffarant Initial soil aoiaturaa, with aaaa albado. 

(6) Ualkar and Rowntraa (1977). 

(7) Rowntraa and Bolton (1978). 

B. Dlffarant albadoa, with aaaa Initial soil molatura. 

(8) Charnay, Quirk, Chow and Kornflald (197/). 

(9) Charvln (1979). 

III. Hybrid Exparlwanta: 

Non-lntaractlva va. Intaractlva soil aoiaturaa. 

(10) Manaba (1973). 

(11) Kurbatkln, Manaba and Hahn (1979). 
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A. Different Soil-Moistures , with Same Albedo. 

(1) ShuVU Mints (1981) . 

The experiment of Shukla and Mints (1961) used the general circulation 
model of tha NASA Goddard Sr>aca Flight Cantor, Laboratory for Ataospharle 
Sclcncas . Tha propartlaa of tha GLAS aodal and Its ability to slaulats tha 
regional and seasonal characteristics of tha observed climate of the earth 
have been described by Shukla at el. (1981). In the experiment, one climate 
simulation is made In which the land-surface evapotranoplratlon, E, la every- 
where made equal to the model calculated potential evapotranspltetlon, Ep, 
which makes the evapotransplratlon coefficient, 6 “ E/Ep - l. \n the other 
case, no land-surface evapotranoplratlon Is allowed to take place at all 
(6 1 0). The prescribed albedo la the same In both cases, and la a very 
slightly modified version of the one given by Posey and Clapp (1964). For 
convenience, the two calculations are called the "wet-soil" case <<nd the 
"dry-soil" case. Both calculations were started from the seme Initial 
observed atmospheric state on 15 June. The results that are shown here 
are the averages for July. 

In the wet-soil caaa, the calculated land-surface evapotransplratlon is 
relatively constant (within about ± 1 mm/ day) between latitudes 20*6 and 60*Mf 
with an average valut of 4.3 mm/ day; corresponding to an evaporative cooling 
of tha surface of 125 watt/m 2 , aa shown in Fig. 2. Kara, tha sensible heat 
u.».ier to tha atmosphere is 21 watt/m 2 . In the dry-soil case, however, the 
lend-eurfaee evapotransplratlon la zero and the sensible heat transfer is 
169 watt/m 2 . 

The dry-soil case gives rlst' to much less cloudiness over the continents 
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R u « (419 + - 393 + ) = 26 R t * (550 ♦ - 461+) - 89 

R N a 172 - 26 = 146 R N = 258 - 89 = 169 


Fig. 2 Surface energy transfers (watt/m*) averaged for Che continents 
between 20*S and 60*N, in experiment of Shukla and Hintz (1981). 

Rs: solar radlational heating of the ground* 

longwave radlational cooling of the ground 
(difference between radiation emitted by 
ground and radiation absorbed by ground. 

R|f ■ (R S “Rl): net (all-wavlenglh) radlational heating of 

the ground. 

LE: latent heat transfer from ground to atmosphere 

(evaporative cooling, of the ground). 

H : conductive-convective heat transfer from ground 

to atmosphere. 
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than the wet-soil case and, as a result, a larger aaount of solar radiation 
reaches and Is absorbed by the ground, 2S8 Instead of 172 watt/a 2 . The 
Increased solar heating of the ground, as well as the elimination of the 
evaporative cooling, aakes the ground varaer; and the higher ground 
teaperature produces a greater long wave radiation ealsslon froa the ground, 
550 Instead of 419 watt/a 2 . The ataosphere also becoaes varaer In the 
dry-soli case, and there Is an Increase In the atmospheric long wave “back 
radiation** to the ground; but because of the reduction In the cloudiness 
the increase in tae back radiation, froa 393 to 461 watt/a 2 , is only 
about half as large as the Increase in the radiation eaitted by the ground. 

The end result of all these large, but partially coapensatlng changes in 
the radiation transfers. Is that there Is only a relatively small change 
In the net (all-wavelength) radlational heating of the land-surface: an 

Increase of only 23 watt/a 2 froa the wet-soil to the dry-soil case. 

The top panel of Fig. 3 shows the global precipitation distribution in 
the wet-soil case. Over most of North America and most of Eurasia the pre- 
cipitation Is with . a abouc 1 aa/day of the local evapotransplratlon. Only 
over southeast China does the precipitation exceed evapotransplraton by as 
much as 4 aa/day. Over South Aaerlca there Is heavy rain near the equator, 
which Is about 2 aa/day greater than the land-surface evapotransplratlon. 

Across Africa, at about 10*N, there Is a band of rain which la about 4 aa/day 
greater than the local evapotransplratlon. On the other hand , across Africa 
at about 25*N, and across Africa and South Aaerlca at about 15*S, the precipi- 
tation la 2 to 3 aa/day smaller than the evapotransplratlon. Thus, although in 
the wet-soil case there Is a fairly uniform transfer of water vapor to the air 
by the land-surface evapotransplratlon, within the tropics and subtropics there 
are convergences and divergences of the water vapor transports by the large- 
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scale atmospheric circulation, which enhance or diminish the precipitation by 
substantial amounts. 

The dry-soil case, shown in the bottom panel of Fig. 3, produces almost 
no precipitation at all over Europe and most of Asia; and over most of 
North America the precipitation is only about a quarter to a half of that 
of the wet-soil case. Over the equatorial part of South America, on the other 
hand, the rainfall in the dry-soil case is about the same as in the wet-soil 
case: i.e., about 6 mm/day; but, now, the water vapor which produces that 
precipitation comes only from the ocean. 

Across north Africa, the ralnband is about 400 km (one model grid 
interval) farther north in the dry-soil case than in the wet-soil case, and 
weaker by 3 to 4 mm/day. The precipitation in the dry-soil case is about the 
same as the amount by which the precipitation exceeded evapotransplration in 
the wet-soil case: which is to say that the convergence of the water vapor 

transport by the atmospheric circulation is about the same in the two cases. 

Perhaps the most surprising difference of all, when comparing the dry- 
soil case with the wet-soil case, is the southward and westward displacement 
of the region of maximum precipitation in southeast Asia. Over Bangladesh, 
the convergence in the water vapor transport from the ocean in the dry-soil 
case more than compensates for the absence of surface evapotransplration. 

It is in the dry-soil case that the calculated precipitation most closely 
resembles the observed summer rainfall of southeast Asia. 

Fig. 4 shows the ground surface temperature. In the dry-soil case, in 
which there is no evaporative cooling of the ground and, because of the 
reduced cloudiness, more solar radiation is absorbed by the ground, the sur- 
face temperatures north of latitude 20°S are about 15° to 30*C warmer than 
in the wet-soil case. 
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As shown In Fig. 2, Che total non-radlat tonal heat transfer to the 
atmosphere (H + LE) Is not greatly different In the two cases (146 vs. 

169 watt/m^); but In the dry-soil case all of this Is sensible heat transfer, 
which Is confined to the planetary boundary layer. In the wet-soil case, by 
contrast, the larger part of the transfer Is In the form of latent heat which 
warms the free atmosphere and not the boundary layer; whether Immediately 
and locally realized by convective condensation and precipitation, or 
realized at some later time and distant place. Thus, there Is a different 
vertical distribution, and sometimes a different horizontal distribution, 
of the dlabatlc heating. This can produce significant differences In the 
thermally forced atmospheric circulation and, by the geostrophlc adjustment 
process, corresponding differences In the horizontal pressure distribution. 

Fig. 5 shows the surface pressure fields. The top and center panels 
show the surface pressures reduced to sea level. In the two cases. The 
bottom panel shows the difference between the two surface pressures, without 
reduction to sea level. It Is here that we see the change In the surface 
geostrophlc wind. Over most of the land the surface pressures are about S to 
15 mb lower In the dry-soil case, which means enhanced cyclonic circulations 
over the continents. 

In the wet-soll case the trough of low pressure across Africa coincides 
with the lntertroplcal ralnband, as may be seen by comparing the top panels 
of Figs. 3 and 5. This Is the same relationship that we see over the tropical 
oceans. But In the dry-soil case the trough of low pressure Is about 400 to 
800 km north of the ralnband; which Is about the same relationship that Is 
found over north Africa, In nature. 

When the surface pressure Is lower over the continents. It must be higher 


over the oceans. Most of the Increase Is In the mid-latitudes of the central 



and western North and South Pacific Oceans. Examination of the vertical motion 


field (not reproduced here) shows that In these ocean regions there la an 
Increased subsiding motion In the dry-soil case. As Figs. 3 and 5 show, 
not only does the Increased sinking motion suppress the oceanic precipita- 
tion, but the accompanying low-level horizontal velocity divergence, by 
generating antlcyclonlc vortlclty, Increases the antlcyclonlc circulation 
In these regions; and geostrophlc adjustment produces the corresponding rise 
In surface pressure. 

Over the Atlantic Ocean, the vertical motion field In the wet-soll case 
shows a band of rising motion and low-level velocity convergence, which coin- 
cides with the band of oceanic precipitation just north of the equator (top 
panel of Fig. 3). But In the dry-soil case there Is sinking motion over all 
of the tropical Atlantic; and there Is no oceanic ralnband near the equator 
at all. Over the eastern half of the tropical Pacific, the same kind of change 
takes place, but It Is not as pronounced. Thus, the change In the land- 
surface boundary condition also produces large changes In the circulation and 
rainfall over the oceans. 
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(2) Suarez and Arakawa (peraonal communication). 

AC Che reviewer’s suggestion, che sane wet-soli (8-0) versus dry-soil 
(8-1) sensltlvlcy experiment wss made with Che UCLA general circulation model 
by M. Suarez and A. Arakawa ( personal communication). (For a 
description of the model, see Arakawa and Lamb, 1977; Arakawa and Suarez, 

1983; and Suarez et al. , 1983.) There are substantial differences between 
the UCLA and the GLAS models, of which the moat Important, Insofar as the 
present sensitivity experiment la concerned, may be the way In which the 
planetary boundary layer and the cumulus convection are parameterized. 

In the UCLA experiment the Integrations for the two cases were started on 
the first day of July, with the Initial state of the atmosphere taken from a 
previous general circulation simulation. The results that are shown here are 
for the 31-day period starting on 16 July. Again, the prescribed surface 
albedo follows Posey and Clapp (1964). 

In the wet-soll case, the calculated land-surface evapotransplratlon was 
relatively constant (within about i 1 mm/day) between 20*S and 60*N, with an 
average value of about 6 mm/day. This Is about 1.7 mm/day larger than In the 
wet-soll case of the GLAS experiment, and la probably a consequence of the fact 
that the UCLA model produces less cloud cover than does the scheme used In 
the GLAS model and, thereby, a greater net radlatlonal heating of the ground. 

Fig. 6 shows that the precipitation In the wet-soll case la about 
611 mm/day, over almost all of extratroplcal North America and Eurasia, and, 
therefore, does not differ from the local evapotransplratlon by more than about 
1 mm/day. Within the tropica, however, the precipitation exceeds the local 
evapotransplratlon by about 10 mm/day over the Indochina peninsula, by about 
3 to 6 mm/day over a few small land areas that are close to the sea (Guatemala, 
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Fig. 6 Precipitation (mn/day) In the wet-soli case (top) and dry-soli 
case (bottom), In the experiment of Suarez and Arakawa (personal communication)) 
(Precipitation greater than 2 mm/day shaded). 
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southern India, southeast China, Columbia, Venezuela and northeast Brazil), 
and by a few mm/day over a large area adjacent to the Somali coast of 
north Africa. These are regions, therefore, of substantial water vapor 
transport convergences. 

In the dry-soil case, shown in the bottom pan»l of Fig* 6, there is 
almost no continental precipitation at all. Only in an east-west band across 
north Africa is there a significant amount of precipitation, 2 to 5 mm/day, 
produced by a convergence ol the water vapor transported from the oceans. The 
axis of this ralnband, at 10*N, is about 1000 km south of the axis of the low 
pressure trough which, in the dry-soil case, is at about 20*N. 

The change in the precipitation over the oceans is very large near some 
of the tropical and subtropical coastlines, and especially where there are 
embayments. In the wet-soil case there are pronounced minima over the Gulf of 
Mexico and the Bay of Bengal (accompanied by pronounced maxima over the adja- 
cent land areas). But in the dry-soil case, the minima are replaced by maxima 
over the ocean embayments. Similarly, along the coasts of Central America and 
northeast Brazil, the land precipitation decreases and the nearshore ocean 
precipitation Increases in going from the wet-soil to the dry-soil case. 

The striking difference between the experiments with the UCLA model and 
the GLAS model is that, except for the Sahel region of Africa, the UCLA 
model produces almost no continental precipitation in the dry-soil case. 

Examination of the water vapor transport field by the investigators showed 
that there are regions, such as northeast Brazil, where within the planetary 
boundary layer there is a large convergence of the water vapor transported from 
the ocean, in the dry-soil case, but no rain. The interpretation they made 
(Suarez and Arakawa, personal communication) is that with dry-soil there is a 
very large diurnal variation of the ground surface temperature, which produces 
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a very large diurnal variation in the depth of the model's planetary boundary 
layer, growing in thickness during the day and collapsing at sunset; and that 
it is this diurnal oscillation which, without producing clouds, transfers the 
water vapor from the boundary layer to the free atmosphere, where the transport 
is divergent. This transfer of water vapor from the boundary layer to the free 
atmosphere by diurnal "boundary layer/free atmosphere mixing” Is not unlike the 
transfer by detralnment from fair-weather, non-preclpltatlng cumulus clouds, 
described in the introduction. The same condition of upper troposphere velo- 
city convergence and middle troposphere subsidence must be satisfied. Here it 
is the result of upper level outflow from the region of intense convective pre- 
cipitation over the adjacent ocean. 
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(3) Mlyakoda and Strlckler (1981) 
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general 


Mlyakoda and Strlckler (1981) used an early version 
circulation model of the NOAA/ Princeton Geophysical Fluid Dynamics Laboratory 
(Mlyakoda et al . , 1969) to make and compare two different sets of 14 day 
numerical weather predlctlona for the northern hemisphere, In July, when diffe- 
rent distributions of the soil moisture availability, 0, were prescribed. 

The surface albedo was fixed and followed Posey and Clapp (1964). The 
clouds were cllmatologlcally prescribed as a function of latitude and height. 

The convective-adjustment scheme was used, In which the moist convective 
heating of the atmosphere and the convective precipitation depend only on the 
relative humidity and on the temperature difference between adjacent levels In 
the vertical. Thus there Is no penetrative convection and, consequently, the 
sensitivity of the convective heating and precipitation to the amount of water 
vapor In the model's planetary boundary layer la not as great as In the cumulus 
convection parameterization schemes of the CLAS and the UCLA models. 

In one case 0 was everywhere set equal to 0.5. Tor the other case, the 
authors sought a more realistic field of 0; and for this they took the ob- 
served normal distribution of precipitation for the antecedent six mctth 
period, February through July, and relabeled the lsohyets as lines of constant 
0, according to the arbitrary function shown In the top panel of Fig. 7 
(Mlyakoda et al. , 1979). No account was taken of the antecedent evapotransplra- 
tlon. Consequently, as the bottom panel of the figure shows, over the northern 
forests and the wet tundra regions of Canada and Siberia 0 was made as low as 
In the subtropical deserts (although. In nature. It Is near the maximum value 
of 1.0). 

Figures 8 and 9 show the differences In the evapotransplratlon, preclplta- 
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Flf. 7 Soil aolature availability, 0 
6 aonth precipitation (Mlyakoda, at al.. 


In July, 

1979). 
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6 - 6(X,$) nlnus the case of 8 i 0.5, when the ensemble average Is taken of 
the last twelve days of three sets of 14-day forecasts. 

Ue see, In Fig. 8, that where 8 Is reduced there Is, In general, a 
reduction In evapotransplratlon and precipitation. The largest reduction In 
evapotransplratlon, of more than 7 mm/day, is over the central part of the 
north African and Aslan deserts; with the axis of the maximum evapotransplra- 
tlon reduction at about latitude 20*N across Africa. But the axis of the 
largest reduction in precipitation Is at about 12*N across Africa (the Sahel), 
where the precipitation decreases by about 12 mm/day. As shown In Mlyakoda 
and Strlckler (1981, Fig. 7), the ralnband of the lntertroplcal convergence 
zone across north Africa does not change Its position, but its magnitude goes 
down from about 20 to 8 mm/day when the land-surface moisture source in the 
Sahara Is eliminated. There Is, obviously, a large change that takes place in 
the water vapor transport. 

In the case where the average land-surface evapotransplratlon became 
smaller, the average ocean evaporation become larger. In spite of that, 
the average ocean precipitation decreased, In agreement with the experiments 
made with the GLAS and UCLA models. Presumably It Is, again, an enhancement 
of the sinking motion over the oceans that suppresses the ocean precipitation. 

If we compare the top panel of Fig. 9 with the top panel of Fig. 8, we see 
that there Is a negative correlation between the change In the land-surface 
temperature and the change In evapotransplratlon. 

When we compare the upper and lower panels of Fig. 9, we see that there Is 
a negative correlation between the change in surface temperature and the change 
In the height of the 1000 mb surface. This Is true even where the land-surface 
Is not at a high elevation. 

This experiment shows cwo Important things: 1) That It does not require 
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an extreme change In the soil moisture availability (0), such as a change from 
1 to 0, In order to produce large changes in the precipitation, temperature and 
motion fields. And, 2) that the Influence of the land-surface evapotransplra- 
tlon on the atmosphere operates very quickly. Mlyakoda (personal communication) 
reports that sizeable differences In the surface temperature and precipitation 
appeared within a few days. 

[In an earlier study of the role of the surface transfers of sensible 
and latent heat In numerical weather prediction, Gadd and Keers (1970, Figs. 

4, 5, 6) showed that even In a very short range (18-hour) prediction for north- 
western Europe and the British Isles, in August, the inclusion of evaporation 
and sensible heat transfer from the land and the sea surfaces made a 

noticeable Improvement in the predicted rainfall over the land]. 
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B. Different Albedos, with Same Soil Moisture. 

ORIGINAL PAGE IS 
OF POOR QUALITY 

(4) Charney, Quirk, Chow and Kornfield (1977). 

This experiment was made by Charney et al. (1977) with the NASA Goddard 
Institute for Space Studies general circulation model, the properties and 
performance of which have been described by Somerville et al. (1974) and 
Stone et al. (1977). 

The three runs that are shown here use the prescribed field of non- 
lnteractlve soil moisture availability, 0, from Stone et al. (1977), who 

assumed that 0 - 2 x (RH - 15)/85, 0 max - 1, where RH is the observed 

normal monthly mean relative humidity of the surface air. The observed 
relative humidities, for July, were taken from the tabulation by Schuts 
and Gates (1972). In the above formulation, 0 - 1 when the relative 
humidity is equal to or greater than 57.5 X. Consequently, 0 was made 

equal to or dose to 1, and the evapotransplratlon therefore equal to or 

dose to the potential evapotransplratlon over most of the land surface 
of the earth. Only in a small region in the western United States and 
across the central Sahara was the prescribed July 0 smaller than 0.5. 

In the run designated as case **28", the ice-free and snow-free continents 
were assigned a surface albedo of 0.14; except that a higher albedo, 0.35, 
was assigned to the regions of the observed northern hemisphere deserts 
(see Fig. 19, below). 

In a comparison run, "3a” , the albedo was changed from 0.14 to 0.35 
in three additional regions, the "Sahel*', "Rajputana", and "Western Great 
Plains", which are adjacent to deserts (Fig. 19). Otherwise everything 
was the same as in case 2a. 

In a separate comparison run, case "4", everything was again the same 
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Fig. 10 


Table I 


as in case 2a, except that the change of albedo, from 0.14 to 0.35, was 
made in three regions that are within the observed rainy and vegetation 
covered areas of the earth. The locations of these regions, called "Central 
Africa", "Bangladesh", and "Mississippi Valley" are given in the table. 

The first numbered column of Table I shows the prescribed soil mois- 
ture availabilities for the six regions. In the Sahel region B - 0.51, and 
in the other regions 0 is 0.78 or more. 

Fig. 10 shows the evapotranspiration and precipitation in the northern 
hemisphere, for case 3a. Over most of the continents the evapotranspiration 
and precipitation do not differ by more than about ± 1 mm/ day. It is only 
over southeast China and Indochina, and where the inter tropical ralnband 
crosses Africa, that the precipitation exceeds the evapotranspiration by 
3 to 4 ram/ day. Near the Mediterranean coast of Africa and in the Middle East 
the precipitation is less than the evapotranspiration by about 2 to 4 mm/day. 

Table I is a rearrangement of the data in Charney et al. (1377, Tables 
4.1 to 4.4) and shows the components of the energy and water budgets at the 
earth's surface in the three desert-margin regions, (case 2a vs. 3a) and in 
the three humid regions (case 2a vs. 4). 

Ue see, in column (10), that in the Western Great Plains (where 0 was 
assigned the value of 0.78) and in the Mississippi Valley (where 0 was made 
1.0), the evapotranaplrations with albedo of 0.14, are, respectively, 4.2 
mm/day (130 mm/mo) and 5.1 mm/day (158 mm/mo). These are in fair agreement 
with the aerologlcally derived evapotranspiration over the central and 
eastern United States in July, of 3.7 mm/day (116 mm/mo), shown in Fig. 1. 
More Important, however, as an indication of the reliability of the model, 
is the fact that the vertical Integral of the water vapor flux convergence. 
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Table I. Components of the energy and water budgets, In experiment of Charney et al. (Ifr77) 
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0 : average, Cor Che region, of Che prescribed soil moisture 

availability (ratio of evapocransplratlon Co poCenClsl 
evapotranspl radon) 

(1 - a) : fraction of Che incident solar radladon Chat is absorbed 

by Che ground (a ■ land-surface albedo) 

2 

R§ : solar radladonal heating of Che ground (vaCC/m ) 

Rl : longwave (Infrared) radladonal cooling of che ground 

(difference becveen longwave radladon emltced by Che 
ground and atmospheric "back radladon" absorbed by Che 
ground) (watt/nr) 


Rtf ■ (R§ - Rl) : net (all-wavelength) radladonal heating of Che ground 

(watt/m 2 ) 

LE : latent heac transfer from ground to atmosphere (evaporative 

cooling of the ground) 

H : conductive-convective heat transfer from ground to atmos- 

phere 


N 


fraction of the sky covered by clouds of all types 


Ta : surface air temperature (*C) 

E : surface evapocransplratlon (mm/day) 


(-V*q\v) : vertically Integrated convergence of the water vapor trans- 

port (ma/day) 


P : precipitation (mm/day) 

AP/AE. : ratio of precipitation change to evapotranspl ration change 


For each region the third line shows either the absolute change between the 
two cases or the percentage change, where the Z sign indicates the latter* 
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Fig. 11 
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(-V«q w), shown in column (11), is negative In the two regions, 
with values respectively of -0.5 mm/ day (-16 ma/mo) and -0.7 am/day 
(-22 mm/mo). This means that water vapor Is being exported from these re- 
gions at about the same rate as the observed transport divergence, of 
0.7 ma/day (22 ma/ao), shown in Fig. 1. 

In the other four regions, in the case with normal surface albedo, 2a, 
the vertical Integrals of the water vapor transport convergence are positive, 
water vapor Is being Imported (so that precipitation Is larger than evapotran- 
splratlon), which Is what one would expect for these particular regions In 
the month of July. 

Tsble I Is replete with Information about the performance of the model 
and its complex, non-linear response to the change In the surface albedo. But, 
for brevity, we will here examine only what happens In the Sahel, the region 
of greatest Interest. 

We see, In Table I (column 3), and In Fig. 11, that when the surface 
albedo Is Increased In the Sahel, from 0.16 to 0.35, the solar radlatlonal 
heating of the surface does not bee one smaller: It becomes larger . This Is 
because of the large decrease In the cloud cover, from 0.70 to 0.66, (column 
8), which more than compensates for the Increased albedo. 

The cloud cover Is less because 1) there Is less evapotransplratlon 
(a change from 3.7 to 2.8 am/day), and 2) there Is less convergence In the 
water vapor transport (a change from 3.7 to 1.2 am/day). 

The local evapotransplratlon Is reduced In the hlgh-albedo case because 
there Is aore long wave radlatlonal cooling of the ground (an Increase from 
58 to 86 vatt/a^). Unfortunately, no record was kept of the ground surface 
temperature, nor of the long wave emission by the ground; but It Is most 
I'kely thst It is the decrease In the downwsrd longwave “back radiation** from 
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R L = (R L + - R l *) = 58 R l = (R L *- R l *> = 84 



Fig. 11 The Sahel region energy budgets (vatt/m 2 ) (cop) and water 
budgets (mm/day) (bottom) , In the experiment with prescribed soli moisture 
availability of Charney et al. (1977). Case 2a Is on the left. Case 3a Is 
on the right. (For definitions of symbols, see legend of Table 1). 







ORIGINAL PAGE IS 
OF POOR QUALITY 


Che atmosphere, as a consequence of Che decreased cloudiness, which 

Increased Che long wave cooling of Che ground. 

The & .nslble heat transfer from Che ground Co Che atmosphere is small. 

In both cases. 

By using arrows Co Indicate when a change In one parameter produces a 
change In another, we can describe what happens In Che Sahel, In this 
experiment, as cbe coupling between a sequence of processes operating 

locally and a sequence ChaC Involves Che large scale atmospheric clrcula- 

(changes) . (changes) v 

Cion. The local sequence Is: albedo — — r radiation r 

(change*) 

evapotransplratlon • ■ > * precipitation. The larger scale sequence Is: 

(changes) „ .... (change*) 

preclpltatlon/condensaclon-heatlng ~ — r large scale circulation * 

(change*) v 

water vapor transport convergence x precipitation. 

The second of these two sequences Is similar (but not Identical) to the 
one In the Charney (1975) hypothesis on the dynamics of deserts. In that 

(change*) . (change*) k 

hypothesis the sequence Is: albedo 1 ■ > surface temperature 7 

(change*) . (change*) v 

large scale circulation ~~ — -• — > water vapor transport convergence > 

precipitation. Here,. the change In evapotransplratlon plays no role. (A 
direct examination of the Charney (1975) hypothesis with a general circulation 
model would consist of a comparison of two runs, In both of which no evapo - 
transpiration Is allowed In the. region of Interest, but is allowed elsewhere, 
and the albedo In the region of Interest Is changed). (See experiment (8), 


below. ) 

The relative Importance of the two sequences of processes, when surfae 
evapotransplratlon does take place, may be seen In columns (13) and (14) of 
Table I. We see. In column (13), that In the Sahel, Rajputana and Central 
Africa, the reduction In the water vapor flux convergence Is between 2.8 to 
4.2 times larger than the reduction In the evapotransplratlon. But In the 


\ 


I 
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Western Great Plains the reduction In the water vapor flux convergence is 
only half as large as the reduction in evapotransplratlon. In Bangladesh 
and the Mississippi Valley things go the other way: increasing the surface 

albedo again decreases the evapotransplratlon; but It Increases the water 
vapor flux convergence. 

As indicated earlier, the response of the large scale precipitation, 
temperature and notion fields to a change In the surface boundary conditions 
(whether soil noisture availability or albedo) will depend on aany factors. 
Of particular lnportance is the horizontal scale and the latitude of tho 
region in which the boundary condition Is changed. Through the geostrophlc 
adjustnent process, the horizontal scale and the latitude determine whether 
the circulation change will be in the vertical plane (snail scale or low 
latitude) or in the horizontal plane (large scale and high latitude). 

When the circulation change is in the vertical plane there la a positive 
feedback on the condensation hea. g , through water vapor transport conver- 
gence. But when the circulation change is in the horizontal plana there 
is a negative feedback on the condensation heating, because then the 
transport removes water vspor, as well as sensible heat, froa the region 
of the condensation heating. 
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(5) Canon and Sangater (1981) 
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Fig. 12 


fable II 


Another experiment in thla category vaa eade by Caraon and Sangater 
(1981) with a low reaolutlon (N20) veralon of the Britlah Meteorological Office 
5-layer general circulation model (Corby et al. , 1977). In both runs, evepo- 
tranaplratlon waa made equal to the calculated potential evapotransplratlon, 

(0 = 1). In one caae the albedo of the anow-free land waa everywhere aet equal 
to 0.1. In the other case, It waa everywhere set equal to 0.3. The regaining 
lower boundary condltlona (sea surface temperatures, sea Ice, and land anow 
cover), aa well as the cllaatologlcally determined longwave radlatlonal heating 
rates, tiers the observed July values. 

Figure shows the precipitation averaged over 90 days (Days 21 to 110 
of Integration), where the top panel la the low albedo case, the center panel 
Is the hlgh-albedo case, and the bottom panel Is the difference between the 
two. 

We see that the high albedo caae has less rainfall over moat of the 
continental areaa: but that over the oceans the rainfall la increaaed. 

The averages of the land precipitation,, and of other parameters, are 
shown In Table II. 

Like the experiment by Charney et al. , (A) above, thla Is also an albedo 
change experiment with permanently wet-soil (0 s 1). Here, too, the Increase 
of albedo produced a decrease In evapotransplratlon (-0.9 mm/day) and an even 
larger decrease In precipitation (-1.2 mm/day); but now the albedo Is changed 
on the continental scale, whereas In (4) It was changed on a acale of only a 
few hundred kilometers In the widths of the various regions. It la not 
surprising, therefore, that In this experiment the contribution to the change 
In precipitation of the change In the water vapor transport convergence la 
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Fig. 12 Precipitation In the low a 
albedo case (center), of the experiment b 
The contours sre for 1, 2, 5, 10 and 20 
< 1 mm/day; heavy shading, pteclp > 5 turn/ 
Is the difference In the precipitation: 
high albedo case (unshaded area positive) 
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only about a third as large as Is the contribution by the change In the local 
evapotransplratlon. Both experiments have about the sane average water vapor 
transport convergence, 0.8 nun/day vs. 0.85 mm/day. But In (4), where the albedo 
was changed from 0.14 to 0.35 over a number of small regions, the average change 
In the transport convergence In those regions was -1.2 mm/day. In (5), where 
the albedo was changed from 0.1 to 0.3 over all of the land, the change In the 
transport convergence was only -0.3 mm/day; showing, again, that the larger 
the horizontal scale the smaller Is the role of the water vapor transport con- 
vergence In compensating for a decrease In the land-surface evapotransplratlon. 


Table II. Components of the energy and water budgets. In the experiment 

of Carson and Sangster (1981 )~ 

Global Averages Over Land, 90— day means (Days 21—110) of a permanent July. 


Sfc. Albedo 


(2) 

H 

(3) 

E 

(4) 

(-7*q \v) 

(5) 

P 

(6) 

A(-V»q \v) 

(7) 

AP 

AE 

AE 

0.1 

104 

35 

3.6 

1.0 

4.6 



0.3 

78 

21 

2.7 

0.7 

3.4 



Difference 

-26 

-14 

m 

B 

B 

D 

m 


For definition of symbols, see Tsble I. 
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II. EXPERIMENTS WITH INTERACTIVE SOIL MOISTURE 
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In all but one of the experiments that follow, the time-dependent soil 
•. moisture is governed by the equations: 


- P - E, W * W* , 

(1) 

' max 


8 E p , 

(2) 

JL , 8 - 1 , 

kW* max 

(3) 


where W is the available moisture stored in the soil, W* is the available 
moisture storage capacity of the soil, P is the rate of precipitation, E is 
the rate of evapotranspiration, Ep is the rate of potential evapotranspira- 
tion, 6 is the soil moisture availability, and k is a prescribed coefficient 
(cf. Carson, 1981). In all of the models, E p is evaporation calculated by 
an aerodynamic method, under the assumption that the vapor pressure at the 
surface is the saturation value for the calculated ground temperature. 


I 
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A. Different Initial Soil Moistures, with Same Albedo 


(6) Walker and Rowntree (1977). 
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OF POOR QUALITY 


Walker and Rowntree (197/) examined the Interaction between time-dependent 
soil moisture and the calculated precipitation, temperature and circulation of 
the atmosphere, not in the global don tin, but in a zonal channel between lati- 
tudes 16 # S and 36°N, and extending over 32" of longitude with cyclic east-west 
boundary conditions. The la;.d and s^a distribution was made zonally-symmetrlc, 
with land to the north and ocean to the south of 6 a N latitude; this being an 
idealization of the western part of north Africa and the Gulf of Guinea. 

The model was an ll-lay.. primitive equations model with 2° latitude- 
longitude resolution. The radiational part of the thermal forcing was taken 
as a constant radiational cooling of the atmosphere, of 1.2*K per day from 
the surface to the 200 mb level, with radiative equilibrium at higher levels 
(which means a constant radiational cooling of the atmosphere of 110 watt/m^); 
and with a constant net radiational heating of the land-surface, R(y, of 150 
watt/m^. Thus, over the land, there was a prescribed horlzontally-unlform 
radiational heating (of AO watt/m?) of the atmosphere-earth system; but 
over the ocean, the surface temperature, and not the surface radiation flux, 
was the prescribed boundary condition. The prescribed, zonally-symmetrlc 
ocean temperatures, from 16°S to 6°N, were the observed August normals at 
0* longitude. The moist-convective adjustment scheme was used to obtain the 
convective precipitation and moist-convective heating of the air. 

The available soil moisture and the land-surface evapotranspiratlon 
were calculated with the equations given at the beglnn'iig of this section; with 
W* taken as 150 ran, and k taken as 0.333. Therefore, 6 ■ 1 when W ^ 50 ran. 
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Two lntegradons were made, in which everything was the same, except 

that : 


In Case I (the initially dry-soil Sahara), W was initialized at zero in 
the latitude zone 14°N - 32°N; and at 100 mra in the land zones 6°N - 14°N and 
32°N - 36°N. 


■ig. 13 



In Case II (the initally moist-soil Sahara), W was initialized at 100 mm 
over all of the land region, 6 C N - 36°N. 

Figures 13, 14 and 13 show the time-evolutions of the soil moisture 
and precipitation in the two cases. 

In Case I (Fig. 13), where the soil was initially dry between 14° and 
32°N, it remains dry. There is almost no net water vapor transport into that 
region and, therefore, there is no precipitation and no water is added to the 
soil there. On the other hand, in the land region 6°N - 14°N the initial 
soil moisture, of 100 mm, goes down to about 90 mm over the first seven days, 
showing an excess of evapotransplratlon over precipitation which averages 
about 1.4 mm/day. The corresponding seven day water vapor transport diver- 
gence, of about 1.4 mm/day, is the difference between a large northward 
transport of water vapor across the coastline by the mean meridional circu- 
lation (Walker and Rowntree, 1977, Fig. 4. a) and an even larger equatorward 
eddy-transport of water vapor by the wave disturbance which developed and 
moved westward across the region. By the end of the Integration period, 
this part of the system also appears to have reached a steady state, except 
for a short period and small amplitude variation produced by transient waves 
in the flow. 

In Case II, the initially moist-soil Sahara, (Figs. 14 and 15), there 
is a rapid development of precipitation in the zone near the coast which, 
after about 2 days, exceeds the evapotransplratlon rate and the soil moisture 
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Precipitation (mm/day) 



Latitude 


Fig. 13 Variation with tine of the tonally-averaged precipitation 
(top) and soil Moisture (bottom) in the case where initially the soil in 
the Sahnra is completely dry. Experiment of Walker and Rowntree (1977). 
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Soil moisture (cm) 


Fig. 14 Variation with time of the zonally-averaged precipitation 
(top) and soil moisture (bottom) In the case where the Initial soil mols~ 
ture In the Sahara Is 100 idb. Experiment of Walker and Rowntree (1977). 
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P 


Fig. 15 Variation of the tonally-averaged precipitation and toll 
moisture with time, In selected latitude zones across north Africa, In 
the case where the Initial soil moisture In the Sahara Is 100 mm. 
Experiment of Walker and Rowntree (1977). 







starts to Increase. The average precipitation in this coastal zone reaches 
18 imo/day on day 5; with a maximum of 30 mm/day at 9*N. After that the 
precipitation rate in this zone decreases rapidly and seems to be starting 
an oscillation about an average rate of around 6 ram/day. The prescribed 
R n - 150 watt/m 2 would provide enough energy for E - E p - 5.2 mm/day; but 
the calculated evapotransplratlon may be smaller or larger than this, de- 
pending on whether the sensible heat transfer at the surface is upward or 
downward . 

Over the rest of the land region, 14° - 32*N, the evapotransplratlon 
exceeds the precipitation until about day 12; and, thereafter, except for an 
oscillation produced by the transient wave disturbances, evapotransplratlon 
equals precipitation and the soil moisture remains constant. 

From what we see in these figures, it appears unlikely that the solutions 
for the initially dry-soil Sahara and the initially moist-soil Sahara will ap- 
proach one another no matter how long the Integrations were to continue. It 
seems safe to say that this highly simplified soil moisture-atmosphere system is 
intransitive. 
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(7) Rowntree and Bolton (1978)» 
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Rowntree and Bolton (1978) made an interactive soil moisture experiment 
with the 5-layer, 500 km grid size, version of the British Meteorological 
Office general circulation model (described by Corby et al. , 1977). 

For the calculation of the soil moisture and evapotransplratlon, 

W* was taken as 200 mm and k as 0.5; so that B ■ 1 when W 2 100 

Three 50 day Integrations were made, all starting from the same initial 
atmospheric conditions on 27 May, but with different initial distributions of 
soil moisture. 

In one run, designated C (for control), the initial soil moisture was 
set at 50 mm at all land points over the globe. Therefore, the initial soil 
moisture availability, 8 , was 0.5 everywhere. 

In the run designated W (for wet-soil case) the initial soil moisture 
was set at 150 mm at all of the European land points that are within the re- 
gion enclosed by the rectangle in Fig. 16; but with an initial value of 50 mm 
at all other land points over the globe. Thus, the lnltal 6 was 1 in the 
European region, but again 0.5 at all other landpolnts over the globe. 

In the run D (the dry-soil case) the soil moisture in the European 
region was initialized at 0 mm (but, again, at 50 n ■ elsewhere). Now, the 
initial 0, and hence the initial evapotransplratlon, was zero in the European 
region. 

Figure 16 shows maps of the three rainfall distributions, averaged for 
the 30-day period 15 June - 15 July; and Fig. 17 shows meridional profiles of 
those time-averaged rainfalls along longitude 13*E, where the maximum rainfall 
occurs. In general the differences between the 30-day rainfalls sre comparable 
to what are produced by the natural variability of a model atmosphere with 
fixed surface boundary conditions. But within the European region the tlme- 
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60W 40W 20W 0 20E 40E 60E 


Fig. 16 Precipitation (aa/day) averaged for 15 June - 15 July, In 
experlaent of Rowntree and Bolcon (1978). Center panel: control run, 

where the Initial aoll aolature, on 27 May, was 50 aa everywhere. Top 
panel: caae where the European land polnta within the indicated 

rectangular region had an initial aoll aolature of 150 aa. Bottom panel: 
caae where the European land polnta within the indicated rectangular 
region had aero aoll aolature. 









averaged precipitation, for daya 20 to 50 following Initialization, was 
greatest in the case which Initially had the most moisture In the soil, 
and smallest when the soil was Initially devoid of moisture. 

The changes In the precipitation and In the soil moisture with time are 
shown In Fig. 18, where the values are 10 day averages for the Indicated 
European land region. Ue see that even after one and a half months there are 
still large differences between the precipitation rates and between the soil 
moistures, when we compare the Initially wet-aoll case (W Initially 100 mm) 
with the other two case (W Initially 50 mm and 0 on). The slopes of the 
rainfall curve and soil moisture curve, for the initially wet-soil case, 
suggest that the system Is transitive; but that the time required for con- 
vergence Is several months. 
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Fig. 16 Variation with time of 10-day averaged precipitation (top) 
and aoll mol a cure (bottom), averaged for the Euopean land polnta within 
the rectangular region ehown In Fig. 16, for the three caaes In the cxperl 
■ent of Rowntree and Bolton (1976). 
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B. Different Albedos, with Same Initial Soli Moisture. 

(8> Charney, Quirk, Chow and Kornfleld (1977). 

Charney et al. (1977) made an experiment in which they compared tiro runs 
which had different time-constant albedo distributions, but the same initial 
soil moistures which could Interact with the atmosphere. 

The distribution of the land-surface albedo in the two cases is shown in 
?lg. 19 the bottom panel of Fig. 19. In their run designated “2b", the albedos of 

“permanent desert" (the regions with dotted shading) were assigned the value of 
0.35; and everywhere else over the ice-free and snow-free land surface of the 
globe the albedo was taken as 0.14. In the comparison run, designated “3b“, 
three regions adjacent to the permanent deserts, the “Western Great Plains", 
“Rajputana” and "Sahel" (shown by the cross-ruled shading) were also assigned 
an albedo of 0.35. In both runs, the initial soil moisture was taken to 
be zero everywhere. 

The change in the time-dependent soil moisture was calculated, in half- 
hourly time steps, from Eqs. (1) and (2) given at the beginning of this section, 
and with the function B - 0 (W,W*,E p ) (Charney et al. , 1977, p. 1368) which la 
?lg. 20 shown in Fig. 20. (Figure by personal communication from Y. Sud.) Over the 
range of E p between 1.4 and 6.4 mm/day, this formulation for 0 was taken from 
Denmead and Shaw (1962), who obtained it from measurements of dally (24 hr) 
evapotransplration and notentlal evapotransplratlon, together with measured 
soil moisture. 

In the experiment, this formula for 0 was Inadvertently applied to the 
calculation of half-hourly values of evapotransplration, and this made the 
calculated dally evapotransplratlon an order of magnitude too small (because 
during the mid-day hours, when E p is about it times as large as its 24 hour 
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. rw Fl8 ‘ 19 BOt 5°“ panel: Che assl 8 ned albedos In the experiments 

all cllet* a nd*dot Unshadcd land areas hav « an albedo of 0.14 i n 

all cases, and dot shaded areas have an albedo of 0.35 In all cases In 

£> rars: ::::: sivsis” changcd fro ° °* u (in 2a - 

Top panel, precipitation (mm/day) In case 3b. 
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average, the 8 obtained In this way is extremely small for almost all values 
of W/W*; and during the night, when 8 can approach 1, Ep Is negligible). 

The top panel In Fig. 19 shows the calculated July precipitation in 
Case 3b. At the beginning of the run, 18 June, almost all of the precipita- 
tion that falls on the land must be from water vapor transported from the 
ocean; and this situation must continue (because evapotransplration is 
negligible) until the accumulation of the water in the soil brings W close 
to the soil storage capacity W*, at which point 8 becomes equal to 1 for 
all values of Ep. For a precipitation of 5 mm/day, which is the mean 
July precipitation across north Africa in Case 3b, it takes about 20 days, 
or until about 8 July, for the soil moisture, W, to approach or reach W*, 
at which point 8 approaches or reaches 1. Elsewhere, 8 remains negligible 
until even the end of July. It is therefore not surprising that in this 
experiment with interactive soil moisture the precipitation, for days 14-44 
following initialization with zero soil moisture, should resemble the per- 
sistently dry soil case shown in the lower panel of Fig. 3. 
ible III Table III shows the energy and water balances of the three regions, in 

the two runs, for July (days 14 to 44 of the integration). We see, by 
comparing columns (11) and (12), that during this period (which, as indicated, 
is a transient stage for the soil moisture and evapotransplration) the domin- 
ant term in the supply of water vapor for precipitation, in all cases, is 
the water vapor transport convergence. In the Sahel, the evapotransplratlons 
for the month are 0.14 and 0.34 mm/day, which are only 4Z and 12Z of the 

monthly averaged precipitation rates. But almost all of this must be due 
to the evapotransplration near the end of July, at which time, or shortly 
thereafter, the accumulated precipitation will make W approach or equal W*. 

With the small amounts of evapotransplration, the cloud cover (column 8) 
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Table III. Components of the energy and water budgets. In experiment of Charney et al. (1977 



ORIGINAL PAfip 

* saw 



















ORIGINAL PAGE IS 

OF POOR QUALITY 

is noc very different in the low and high alhedo cases. Consequently, unlike 
Experiment (A) above, the change in the net radlatlonal heating of the ground, 

Rjq , is almost entirely due to the change in the solar radiational heating 
of the ground, Rg; and this produces the large change in the sensible 
heat transfer, H. It is this change in the sensible heating of the planetary 
boundary layer which produces the change in the vertical velocity shown In 
Fig. 21 the bottom panel of Fig. 21; the associated change in water vapor transport 
convergence, shown in column (11) of Table II; and, the change in precipita- 
tion shown in column (12) and in the top panel of Fig. 21. As the experiment 
stands, it does illustrate the mechanism of the Charney (1975) “dry-soil'* 
desertification hypothesis: but not if the integrations were to continue; 

for then W will everywhere approach or equal W* and the processes which 

depend on evapotransplratlon will become Important. ^ 

We note, furthermore, that this simple picture of tne coupling of the 
water vapor transport convergence to the surface albedo does not hold for the 
Rajputana and Western Great Plains regions, even when the evapotransplratlon 
is negligibly small. In both of these regions the changes in Rg, R^, and H 
are about the same as the changes In the Sahel; but, unlike the Sahel, there 

r 

is almost no change in the water vapor transport convergence. As already 
indicated, the response of the atmospheric circulation to a change In the 
boundary layer heating, H, will depend very much on the horizontal scale and 
latitude of the heating perturbation, as well as on the orientation of the 
region with respect to external moisture sources and the way In which the 
altered flow encounters the mountain barriers. 
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Fig. 21 Zona lly-ave raged precipitation (top) and vertical velocity 
In the alddla troposphere (bottoa), over Africa, when the albedo In the 
Sahel (l6*N - 20*N) Is Increased froa 0.14 co 0.35. Experiment of Charney 
et'-al. (1977). 
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(9) Chervin (1979). 
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Chervin (1979) used the NCAR general circulation model (described by 
Washington and Williamson, 1977) to examine the effect of a change In the 
land-surface albedo when the soil moisture Is fully Interactive. The change 
In soil moisture was calculated with the equations given at the beginning 
of this section, with W* - 150 mm and k - 0.75. 

The control was the average of a master run, which started from a state 
of res*: and isothermalcy and was Integrated for 120 days, plus four other 
runs, each of which started from day 30 of the master and ran until day 120. 
All of these were perpetual July integrations, In which the sun declination, 
the ocean surface temperatures, and the snow-free land albedo (which followed 
Posey and Clapp, 1964) were held constant In time. 

The run with a different albedo was also started from day 30 of the 
master and ran until day 120. The change in the albedo consisted of replacing 
the Posey and Clapp values by a constant albedo of 0.45 within two regions: 

1) a large region over north Africa, extending from the Atlantic to the Red 
Sea and from latitude 7.5' J N to the Mediterranean; and therefore covering the 
zone of the July lntertroplcal convergence rain as well as the Sahara Desert; 
and 2) a smaller region over the U. S. High Plains (97.5 # W - 107. 5*W, 

27 .5*N - 52 . 5*N) . 

Over Africa, the control (Posey and Clapp) albedo varies from about 0.35 
In the northern Sahara to about 0.08 near the southern boundary of the region 
where the albedo will be changed. Over the U. S. High Plains, the control 
albedo Is between 0.07 and 0.17. 

Fig. 22 shows the change over Africa In the precipitation, soil moisture 
and ground temperature, and In the vertical velocity at 3 km elevation. The 
values shown are the averages of the last 60 days of the modified albedo case, 
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A 3KM VERTICAL VELOCITY (MM/SEC) A PRECIPITATION RATE (MM/DAY) 



A GROUNO TEMPERATURE (*C) A SOIL MOISTURE (CM) 


Fig. 22 , The change in verdcal velocity at 3 km elevation (top left), 
precipitation (top right), ground temperature (bottom left), and soil mois- 
ture (bottom right). In the albedo change experiment of Chervln (1979). 
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minus the ensemble average of the last 60 days of the five control runs. The 

stippled areas in the figure show the regions of r - |A |/o > 3; where 

60 60 “ 

A is the prescribed change response (i.e., the difference between the 60-day 
60 

mean in the prescribed change case and the ensemble average of the 60-day means 

of the five control runs); and a is the standard deviation of the 60-day means 

60 

of :he five control runs. According to Chervln and Schneider (1976), r _> 3 

implies an approximately 5X significance level in rejecting the hypothesis 

that the prescribed change response is the result of random fluctuations and 

not the result of the prescribed surface albedo change. 

The maps in Fig. 21 show that the changes are greatest at and near the 

zone where the albedo is Increased from 0.08 to 0.45. In this zone there is 

a decrease in the average upward motion of the air of about 2 mm/sec 

(200 m/day), with r > 3; a decrease in the average precipitation of about 

4 mm/day, with r > 2, (A P - (4 - 8) mm/day, a P - 2 mra/day); and a 
■" 60 60 

decrease in the average soil moisture storage of about 50 mm, with r > 3. 

There is a decrease of the ground surface temperature, with r > 3, over 
almost all of the region of the albedo change, but not along its southern 
edge. There, the ground surface temperature increases , by about 0.5*C, with 
r > 3. 

The paradoxical rise in the ground surface temperature, in the region of 
the largest increase of albedo (the change from 0.08 to 0.45), can be attri- 
buted to the fact that in that zone, 7.5*N - 12.5°N, where there Is the largest 
decrease in precipitation, the soil moisture, U, goes down from about 100 am in 
the control to 50 mm in the high albedo case. With W* ■ 150 mm and k ■ 0.75, 
this reduces 8 from 0.9 to 0.45, and, consequently, there is a large reduction 
in the evapotransplratlon and the evaporative cooling of the surface in that 
zone. There Ip also a contribution to the temperature rise from the accompanying 
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reduction of the cloudiness In that zone (R. Chervln, personal communication). 

Over the U. S. Great Plain region there was almost no change In the 
vertical velocity at 3 km elevation; an average decrease of about 1 mm/day 
In the precipitation, with r > 3 over about half of the region; almost no 
change In the soil moisture storage; and a decrease, averaging about 2*C, 
with r > 3, In the ground surface temperature. 
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III. HYBRID experiments 
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Non-Interactlve vs. Interactive Soli Mol»ture 

These are experiments In which a calculation with Interactive soil moisture 
is compared with one in which the prescribed soil moisture la held fixed for the 
duration of the experiment. If we regard the Interactive case as a simulation 
of nature, then the case with the prescribed, fixed soil moisture can be thought 
of as showing how the climate would be changed if the land-surface evapotranspl- 
ratlon were to be brought under man's control: as, for example, by large scale 

irrigation or by a change or a complete removal of the vegetation cover. 

(10) Manabe (1975) 

A massive irrigation simulation experiment was made by Manabe (1975) with 
one of the GFDL general circulation models. The model used the moist-convective 
adjustment method for calculating the convective precipitation and the moist- 
convective heating of the air; solar and long wave radiation transfers calculated 
with a non-interactive cloud diatributlon, prescribed as a function of latitude 
and elevation; and an albedo for ice-free and snow-free land that follows Posey 
and Clapp (1964). 

In the "natural-case” (the interactive soil moisture case) the soil mois- 
ture was governed by the equations given at the beginning of Section II, with 
U* - ISO mm, and k - 0.75. In the "irrlgation-caae" (the non-lnteractlve soli 
moisture case) 0 was everywhere held equal to 1. 

The "natural-case" simulation produced a ralnband across north Africa in 
in which, averaged between 15*E and 30*E, there was a rainfall maximum of 
about 6 mm/day at latitude 5*N. In the irrigation-case, the maximum rainfall 
«»*s about 12 mm/day at latitude 8*N (see Manabe, 1975, Fig. 3). 
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Ti&i_ 23 The sol id line in Fig. 23 shows the change In the precipitation, the 

Irrigation-case nlnus the natural-case. The dashed line ahows the corresponding 
difference In the evapotransplratlon. (The reduction In evaporation, between 
.* 31 *N and 37*N, Is over the Mediterranean Sea; presuaably because the air was 

more humid from the massive land irrigation). 

Between 18*N and 30*N the Increase in precipitation Is somewhat less than 
the Increase In evapotransplratlon, while between 12"N and 15*N, where the 
evapotransplratlon Increases by only about 1.5 mm/day, the precipitation goes 
up by about 5.5 mm/day. On the other hand, at the equator, where there Is an 
Increase of evapotransplratlon of 0.7 mm/day, the precipitation decreases by 
3.5 mm/day. It Is obvious, therefore, that there are large changes In the 
water vapor transport convergence. 

"lg. 24 Figure 24 shows the circulation In the meridional plane, (v J + w Ik', 

averaged between 15*E and 30*E (where v Is the northward component of the 
horizontal velocity, w Is the vertical velocity, J Is unit horizontal vector 
directed northward, and k Is unit vertical vector directed upward). Although 
the eastward component of the horizontal velocity does not appear here. Its 
divergence, du/3x, enters Into the calculation of the vertical velocity, w. 

Ue see that at the equator there Is a change from upward motion to down- 
ward motion In the free atmosphere; which must be accompanied by c change 
from horizontal velocity convergence to horizontal velocity divergence In 
the boundary layer. It therefore Is the decrease In the boundary layer water 
vapor transport convergence (-V*qv * " qV» v) which makes the precipitation 

e , 

0 

decrease by 3.5 am/day (from 5 to 1.5 mm/day). 

Between 12*N and 15*N, on the other hand, weak ascending motion changes 
to ve.y strong ascending motion; and the accompanying large Increase In the 
boundary layer water vapor transport convergence, added to the small Increase 
In evapotransplratlon, makes the precipitation Increase by 5.5 mm/day (fro* 

0.7 to 6.2 mm/day). 
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Fig. 23 Change In evapotranspl ration (broken line) and In preclplcadon 
(solid line), averaged between 15*E and 30*E, In the hybrid experiment of 
Manabe (1975). 
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Fig* 24 Circulation in the meridional plane, averaged between 15*E and 
30*E ( in the hybrid experiment of M ana be (1975). Top: interactive soil 

moisture case* Bottom: ' moist soil case* (figure by personal communication) 
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(11) Kurbatkln, Manabe and Hahn (1979). 

The Irrigation scheme of Experiment (10) (which makes 3 ■ 1 everywhere) 

• • is physically realizable in the real world. But it is not so easy to make 

» 3-0 everywhere in the real world. If the rainfall were to be very constant 

in time, so that the surface of the earth was always wet, then it would be 
very difficult (although possible) to prevent evaporation. But where rainfall 
occurs in an Intermittent way, and most of the water infiltrates to a depth of 
more than a few centimeters, then the removal of the vegetation would stop the 
transpiration and, thereby, would greatly reduce the transfer of water vapor 
to the air. 

Kurbatkln, Manabe and Hahn (1979) made a hybrid experiment, in which a 
simulation with 0-0 everywhere was compared with a simulation with inter- 
active soil moisture. For this experiment they used the M-21 version of the 
GFDL spectral model (Manabe, et al. , 1979), with moist-convective adjustment, 
prognostic clouds, and a prescribed albedo that follows Posey and Clapp (1964). 
In the interactive soil moisture case, the soil moisture and evapotranspira- 
tion were calculated with the equations given at the beginning of Section II, 
with W* - 150 mm, and k - 0.75. 

The Integration in the Interactive case was over a period of two years 
and eight months. The results that are shown here are averages for the months 
of July and August at the end of the integration period. The non-interactive 
case, with 0-0, was initialized from the interactive case at the beginning 
of the last June and run until the end of August. 
g. 25 The top panel in Fig. 25 shows the simulated precipitation in the 

interactive soil moisture case, the center panel shows the precipitation 
in the no-evapotransplratlon case, and the bottom panel shows the difference 
between the two (no-evapotranspiration case minus interactive case). With 
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•Fig. 25 Precipitation (ca/day) averaged Cor July and August, in the 
• experlnent of Kurbatkln et al. (1979). Top: interactive aoll moisture 

case. Middle: no land-surface evapotranspl ration case. (Contours 0.05, 

0*1* 0*2, 0.5, 1.0, 5.0 cn/day. Dotted shading: predp. < 0.1 ca/day; 

ruled shading: predp. >0.5 ca/day).. Bottoa: no-evapotransplratlon 

case alnus Interactive case (shaded area negative). 
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no land-surface evapotransplration, there Is less precipitation over the 
continents and also less precipitation over most of the oceans; with the largest 
decreases , of up to 5 ram/day, over India, the north Indian Ocean and over the 

* 

western part of the tropical North Pacific. Only over the aid-latitude east 
coast of Asia and the adjacent Pacific, and in some longitudes along the equator 
(but not over the Atlantic Ocean sector), is the precipitation larger in the 
8 - 0 case, and by as much as 2 mm/day. 

rig. 26 Fig* 26 shows the sea level pressure fields. We see that in the no- 

evapotransplratlon case the pressure is lower over most of the continental 
areas. But over north-central Asia, over the northeast Atlantic, and espe- 
cially over the extratropical central and western North and South Pacific 
Oceans , the pressures are higher in the no-evapotranspiration case. Although 
only the 8*0 case has the same evapotransplration boundary condition that 
was used in the experiment with the GLAS model (Experiment 1 above), there 
are many correspondences, as can be seen by comparing the bottom panels 
Fig. 25 and Fig. 6, especially over the central and western North and South 
Pacific and over the South Atlantic. 

The authors show (Kurbatkin et al. , 1979, Fig. 4) that in the 8*0 
case there is, on the average, an Increase in the large scale ascending motion 
in the free atmosphere over the continents; and they remark that in spite of 
this Increased ascending motion there is less cloudiness over the continents; 
and therefore a greater absorption of solar radiation by the ground and a 
higher surface temperature than in the interactive soil moisture case. The 
decrease in the precipitation over the oceans, they point out, is produced 
by the increase in descending motion (or weakening of ascending motion) over 
the oceans. 
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Fig. 26 Sea level pressure (mb), averaged for July and AugusC, In 
Che experiment of Kurbatkln et al. (1979). Top: interactive soil moisture 
case. Middle: no land-'surface evapotranspiratlon case. Bottom: no- 

evapotransplratlon case minus interactive case (shaded area negative). 
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SUMMARY AND CONCLUSIONS 


All of the experiments show that the model simulated climates are sensitive 
to the land-surface boundary conditions which affect evapotransplration. When 
soil moisture availability or surface albedo are changed regionally (or 
globally), changes in the precipitation, the temperature and the motion field 
of the atmosphere take place over the corresponding region (or over the globe), 
which are clearly above the level of the natural variability of the model simu- 
lated climates (those which are caused by the shear-flow Instabilities of the 
atmosphere) . 

[Whether a change in the boundary condition of a given region has a 
significant effect on the climate of some distant region is not known. 

Such atmospheric "teleconnections" are not self-evident in these experiments: 
and, like other kinds of forcing from a distance, Lhey require some kind 
of statistical analysis to separate the signal from the background noise]. 

Not only are the regional (local) Influences large and consistent from 
one experiment to another; they are also easily understood, in physical terms, 
after the analysis has untangled the non-linear interactions. Thus we find, 
for example, that under some circumstances an Increase in the surface albedo 
reduces the cloudiness and, in that way, Increases (not decreases) the ground 
temperature. Or, as another example, cutting off the evapotransplration over 
Asia increases (not decreases) the Indian monsoon rainfall, because it changes 
the orientation, with respect to the Himalayan mountain barrier, of the moist 
boundary layer air stream from the ocean. It is this kind of behavior, easily 
understood after the fact, which is very difficult to anticipate beforehand. 

The magnitudes of the changes in climate which are produced by modifying 
the soil moisture availability or the surface albedo, as shown in these experi- 
ments, are about as large as the changes produced by the seasonal change in 



the declination of the sun (and, the reviewer hazards the guess, larger than 
the changes produced by the observed seasonal changes in ocean surface tempera 
tures and extent of the sea-ice). It is very likely, therefore, that the land 
surface evapotranspiration process, whose time scale depends on the magnitude 
of the soil moisture storage relative to the difference between the precipita- 
tion and evapotranspiration rates, is the most important boundary process that 
can produce anomalies in the time-averaged state of the atmosphere (changes 
in climate) on the monthly, seasonal and annual time scales. 

A shortcoming of all existing general circulation models is that they 
calculate the potential evapotranspiration with a formulation which is 
appropriate for evaporation from an open water bucket, but not from a 
vegetated surface (and, especially, not from tall (forest) vegetation.) 

(See, for example, Shuttleworth and Calder, 1979; and Sellers and Lockwood, 
1981). But, hopefully, this will be corrected and future models will take 
into account the vegetation Influence on the water and energy transfers to 
the atmosphere in a realistic way. 

• ••••••••••••• 

The reviewer thanks the several scientists who provided unpublished or 
supplementary information about their experiments; Dr. J. Shukla and Dr. Y. 

Sud who made helpful comments on the review; and Ms. J. Reckley for her 
extraordinary patience and kindness in the typing, and re-typing, of the 
manuscript . 
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Since this review was written four additional papers have appeared, 
all of which show that the atmosphere Is sensitive to the land surface 
evapotransplratlon . 
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